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Abstract. We calculate the conductance of tubular-shaped nanowires having many potential scatterers at 
random positions. Our approach is based on the scattering matrix formalism and our results analyzed 
within the scaling theory of disordered conductors. When increasing the energy the conductance for a big 
enough number of impurities in the tube manifests a systematic evolution from the localized to the metallic 
regimes. Nevertheless, a conspicuous drop in conductance is predicted whenever a new transverse channel 
is open. Comparison with the semiclassical calculation leading to purely ohmic behavior is made. 

PACS. 73. 20. At Surface states, band structure, electron density of states - 73.21.b Electron states and 
collective excitations in multilayers, quantum wells, mesoscopic, and nanoscale systems 
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1 Introduction 

Electronic transport in nanostructures is often affected by 
the presence of impurities. Understanding impurity in- 
duced disorder has been a major concern in the field of 
quantum transport for many years. A vast literature on 
this problem has accumulated since Anderson's discovery 
of the phenomenon of weak localization [lj and the pro- 
posal of the scaling theory [2]. A lengthy introduction to 
the field is out of the scope of this work but the reader is 
addressed to Refs. [3, 4, 5, 6 J for reviews on the topic. 

In this paper we address a specific nanostructure ge- 
ometry, namely a two-dimensional (2D) electron gas on 
the surface of a tube. Our motivation to study this par- 
ticular nanostructure shape is mainly due to the similar- 
ity with rolled-up semiconductor quantum wells recently 
fabricated [7,8,9, fD|ITT] and studied [12,13,14J. Here, it 
is also worth mentioning the apparent resemblance with 
electrons in carbon nanotubes, although these latter sys- 
tems are generally of smaller dimensions. We also stress 
that this geometry is particularly adapted to the theo- 
retical modeling of random impurities. Indeed, by means 
of longitudinal translations and rotations, i.e., the tube 
symmetries, one can always relate the position of two im- 
purities on the tube surface. 

We have implemented the numerical modeling of the 
disordered nanotubes using the scattering matrix formal- 
ism. In this approach, the solution of the Schrodinger 
equation is obtained by treating each scatterer as a black 
box characterized by its transmission and reflection ampli- 
tudes. The solution to the many impurity problem is then 
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found by adequately composing the single impurity scat- 
tering matrices. Our approach is similar to that of Cahay, 
McLennan and Datta for planar wires [15] . A practical dif- 
ference, however, is that we do not present the problem as 
a matrix recursion but, instead, as a global sparse linear 
system. 

When increasing the number of propagating transverse 
modes the system behavior evolves from purely one di- 
mensional (ID) in the one-mode limit towards 2D in the 
infinite- mode limit. In this paper we focus on the few- 
mode regime, or quasi- ID limit, monitoring the evolution 
of the conductance for energies allowing up to five prop- 
agating modes. In all cases the disordered wire is charac- 
terized by the asymptotic exponential localization of the 
wave function. However, it will be shown below that de- 
pending on the energy remarkable variations of the local- 
ization properties are predicted. 

The scaling theory of disordered wires implies that, 
given a wire length L, wire conductance is a random vari- 
able G, characterized by a certain probability distribution 
Pl^{G). Randomness appears due to the infinite ways in 
which the number of impurities Af can be arranged (each 
particular arrangement is usually referred as a disorder 
realization). Besides L, the probability distribution is also 
characterized by parameter 7, usually known as Lyapunov 
exponent [4]. The relation between the Lyapunov expo- 
nent and the mean conductance strongly depends on the 
wire length. In other words, the form of Pc >7 (G) changes 
strongly with L. The precise definition of 7 involves aver- 
aging a function of the random conductance over all dis- 
order realizations in the limit of a very long wire; namely 



7 = lim (7) 
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where we have defined 
1 
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In Eq. (|2j), Go = 2e 2 /h represents the conductance quan- 
tum and M is the number of propagating modes in the 
wire. For a given L and disorder realization the system is 
characterized by a single value of G and its associated 7. 
Changing the disorder realization these variables acquire 
the statistical meaning implied above. Clearly, 7 is then a 
random variable whose mean, in the long wire limit, gives 
the Lyapunov exponent. 

The inverse of the Lyapunov exponent gives the local- 
ization length £ = 7 -1 . This length is physically relevant 
since it specifies how G and 7 are distributed for a given 
L. Two limiting regimes are known. 

a) L ^> £ (localized regime): 7 is normally (Gaussian) 
distributed with mean value 7 and width af oc = 2j/L. 
In this limit we then have 



1 



2lT(J\ 



-(7-7)V2-L^I 

dG 



(3) 



where 7 as well as d^/dG can be obtained from G using 
Eq. ©. 

b) L <C £ (metallic regime): G is normally distributed 



2lT<J B 



with mean value given by 



G 



Go 



e lLl 



(4) 



(5) 



and a constant variance. Indeed, universality in the 
value of cr met is the reason why this limit is also dubbed 
the universal conductance fluctuation regime. In our 
case in which time reversal and spin rotation symme- 
tries are fulfilled this value is a me t = Go [3J . Notice 

that in Eq. (jH) we have considered the possibility that 
the mean value jl = (7) may not be fully converged 
yet to the Lyapunov exponent for a given L <C £. 

Purely Ohmic behavior is characterized by a strict lin- 
earity of the resistance with the wire length. Such a be- 
havior is obtained in a semiclassical description, where any 
possibility of quantum interference is neglected, and the 
scatterers are composed incoherently. Total transmission 
in this case does not depend on the impurity positions 
[T6] , but only on its number J\f which is fixed for a certain 
length L. The semiclassical conductance G SC (L) is then 
fully deterministic and can be written in terms of a single 
parameter a as 



G SC (L) 



Go 



(6) 



Deep in the metallic regime the mean quantum conduc- 
tance also has an Ohmic behavior, as seen immediately 



X 



^0 



di RAN 



a i-\-l 
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Fig. 1. Parameters of the impurity distribution. 



from Eq. ([5j) assuming 7^ = 7 and expanding e lL « 
I + 7L. Notice, however, that this Ohmic dependence may 
be different from the semiclassical one since we may have 
a 7^ 7, as will be shown below for particular cases. Includ- 
ing the next-order contribution from the exponential and 
assuming many propagating modes such that <C jL w 
aL C 1 we find 



<G)-G SC (L) = G 



(aL) (7L) 



(?) 



If a = 7, Eq. predicts (G) = G SC (L) - G /2, but this 
need not be true in general. 

In this paper we calculate the localization length of 
disordered tubes and its evolution with energy covering 
from one to five propagating modes. The above prediction 
of scaling theory in the localized and metallic regimes is 
shown to be fulfilled and the deviation from it in the in- 
termediate regime is explicitly shown. Given a number 
of propagating modes we find a steady increase of local- 
ization length with energy. Quite remarkably, however, a 
large drop is obtained each time a new mode becomes 
propagating. As pointed by Thouless [17] the localization 
length is related to the density of states and, therefore, 
we may expect some type of discontinuity in our quasi ID 
system. This variation in localization length is the cause 
of similar drops in the mean conductance of a disordered 
tube of length L. In this case, increasing the energy the 
tube evolves from a more localized towards a more metal- 
lic regime. 

The paper is organized as follows. In Sec. 2 we present 
the model and our approach to the many impurity scat- 
tering problem (in appendix we also discuss the solution 
of the single impurity problem) . Section 3 contains the re- 
sults and its discussion while Sec. 4 draws the conclusions 
of the work. 



2 Model 

We assume the electrons move on the surface of a cylinder 
of radius p. Disorder is represented by repulsive potential 
barriers mimicking the effect of impurities on the tube sur- 
face. The potential is assumed of a short range (Gaussian) 
type 

V^(r) =Voe~^ 2 ^ 2 , (8) 

where is a specific impurity position. M such impu- 
rities, a rather large number, are present with positions 
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Fig. 2. Sketch of input and output amplitudes of the scatter- 
ing matrix approach for a single impurity (a) and for many 
impurities (b). In the latter case input and output amplitudes 
of successive impurities are related. 

{r^, i = 1, . . . , J\f}. All impurities are assumed to be iden- 
tical, although we would not expect qualitative differences 
in the results presented below if some variation in Vb and 
a was considered. 



2.1 Impurity distribution 

Impurities are assumed to be well separated compared 
with the potential range a. A sequential random distri- 
bution is then represented by the two parameters do and 
d\ (see a sketch in Fig. 1). The first impurity is arbitrarily 
positioned at (z\,0\) = (0,0) and successive ones at 



Zi = Zi-i + d 
0i 2tt RAN 



di RAN 



(z = 2,...,A0, (9) 



where RAN is a standard random number uniformly dis- 
tributed between and 1. We define di = Zi — the 
longitudinal separation between impurities i and i — 1. 
Obviously, on average (di) = do + O.bdi and the total 
length L of the disordered tube is accurately given by 
L = (do + 0.5di)jV. The idea behind this model is that 
all physical effects on electronic motion by an individual 
scatterer have vanished within a distance do around the 
scatterer's position, the electron propagating then freely 
a distance diRAN to the proximity of the next scatterer. 



M 

+ E ^rn(0)e- iMz ~ Z ^ 

171=1 



(10) 



where M is the number of propagating transverse modes 
{4 > m{0)}- In Eq. ([TO]) we have introduced a contact label 
c = r, for left and right asymptotic regions with bound- 
aries at zf* and Zr\ respectively. We have also defined 
S£ = 1 and s r = —1, and the mode wavenumber k rn (see 
the appendix). As it is normally defined [18J, the scat- 
tering matrix relates output current amplitudes to input 
ones 



' mn I- mn 

t W r /(0 




(ii) 



where v m = Tik rn /m* is the velocity of mode rn. 

Assume now the many impurity scenario of the pre- 
ceding subsection. As hinted in Fig. 2b, it is clear that 
input and output amplitudes of successive impurities are 
related. Indeed, 



_ h (i-l) ik m (di-d ) 
a £,m ~ u r,m e J 

Ai) _ h (i+l) Jkm(d i+1 -d ) 



(12) 
(13) 



Notice that the distance between impurities is not a con- 
stant, di = Zi — Zi-i, and that the phases in Eqs. ([T2]) and 
([T3]) are due to the different definition of the asymptotic 
boundaries for each impurity. 

Repeatedly using Eqs. (fTTjh ([T2]) and (p~3|) for the suc- 
cessive impurities we would obtain a linear system of equa- 
tions relating input and output amplitudes of the global 
disordered region, i.e., between first (leftmost) and J\f- 
th (rightmost) impurities. This total transmission corre- 
sponding to, say, unit-amplitude incidence from the left in 
transverse mode mi n to the right in mode m ou t is simply 



t = b w 



(14) 



The linear system yielding the output 6-amplitudes reads 



2.2 Scattering matrix approach 

Let us consider a single impurity, the i-th one, say, and 
write the wave function far from it (outside the region of 
length do? see Fig. 2a) as 

M 
m=l 
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Equation ([T5|) is a highly sparse linear system of 2MM 
equations. It can be very efficiently solved with sparse 
linear solvers like ME48 of the Harwell library [19] . 

As mentioned in the introduction, a major simplifying 
property derived from the tube symmetry is the trans- 
formation relating any two impurities. In fact, defining a 
reference impurity at (zo,#o) the transmission coefficient 
for the i-th impurity, positioned at (^,0$), reads 

/(*) - P -i(o t -e )(\m-Xr l ) f (o) / lfi \ 

where A m is the angular momentum of mode m. The phase 
in Eq. (fT6|) is due to the rotation needed to link the two 
impurity positions. The longitudinal translation does not 
introduce any additional phase here because our defini- 
tion of the scattering matrix implies a translation of the 
reference boundaries for each impurity. This type of phase 
does appear however in Eqs. ([T2]) and (jT3|) . It can be eas- 
ily shown that an identical phase to that in Eq. ([T6]) is 

relating r^ } n , t f( £ n and r /( ^ n with r^, t'$ n and r'^, 
respectively. The overall result is that only the reference 
impurity scattering coefficients are necessary to set up the 
linear system Eq. ([T5]) . 

For a given disorder realization we obtain the linear 
conductance G from Eqs. ([T5]) and ([T4|) using Landauer 
formula to relate conductance with transmission 

G = Go5> m „| 2 . (17) 

mn 

As well known, Eq. (fT7j) corresponds to a two terminal 
measurement and it includes the contact resistance, i.e., 
in absence of any disorder it yields G — MGq. 



3 Results 

Figure [3] displays the results obtained for a tube having 
a fixed number of impurities Af = 1000 and an energy 
range including up to 5 propagating modes. We have fixed 
do = 0.5p and the different symbols are for impurity distri- 
butions characterized by d\ = 0.5p (solid) and d\ = 2.5p 
(open). The mean conductance has a conspicuous saw- 
tooth wave behavior, steadily increasing with energy until 
a new mode becomes propagating and a sharp drop oc- 
curs. The variance, AG = ((G 2 ) — (G) 2 ) 1 / 2 , shows similar 
drops and a clear convergence with increasing energy to- 
wards the universal conductance fluctuation value. Except 
for the first plateau [20] the results for d\ = 0.5p and 2.5p 
look quite similar. 

As mentioned in the introduction, the key to analyze 
the results of Fig. [3] is the localization length t. Figure [H 
displays this variable for the same parameters of Fig. [3l 
We have found £ by disorder averaging the Lyapunov ex- 
ponent in a long wire. More specifically, we have done the 
calculations for 10 4 impurities, for which the wire length is 
much larger than £ for all energies shown in Fig. [H There 
is a qualitative similarity between Figs. [3] and [U The en- 
ergy dependence of the conductance for a fixed wire length 
can be interpreted in terms of the strong variation of the 




12 3 4 



Energy (h 2 /m*p 2 ) 

Fig. 3. Mean conductance (upper) and variance (lower) for 
Af — 1000 impurities. Solid and open symbols correspond 
to di = 0.5p and 2.5p, respectively. Other parameters we 
have used are do = 0.5p and the impurity potential given by 
a = 0.05/9 and Vb = 80 % 2 /(m*p 2 ). The solid line in the up- 
per plot corresponds to the semiclassical model and the dash 
on the right vertical axis of the lower plot gives the universal 
conductance fluctuation value [3]. 




Energy (h 2 /m*p 2 ) 

Fig. 4. Localization length, scaled by the mean impurity- 
impurity distance, for the same parameters of Fig. 3. The solid 
line shows the corresponding semiclassical result for a -1 . 

localization length with the energy. In general, as the en- 
ergy increases the wire evolves from the localized towards 
the metallic regime until the appearance of a new mode 
causes a sudden drop back towards the localized case. 

The first plateau, with only one propagating mode, 
shows a qualitative difference from the others. Remark- 
ably, the mean conductance with 1000 impurities is close 
to the maximal value Go when d\ = 0.5 p. For d\ = 2. hp it 
also shows a different energy dependence with a faster en- 
ergy increase than for other plateaus. This behavior can be 
attributed to the effective one dimensionality in this limit. 
Indeed, it can be shown that the model is equivalent to 
a Kronig-Penney lattice in which the spacing between ID 
barriers is not exactly regular but has some fluctuation 
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around the mean positions. For a fixed do regularity is 
inversely proportional to d\. The more regular the distri- 
bution, the higher the conductance due to the formation 
of extended Bloch states. This is nicely confirmed by the 
localization length (Fig. [4)) which shows a dramatic en- 
hancement in the first plateau for the smaller d\. How- 
ever, for large enough wire lengths the states are always 
localized, i.e., £ is large but remains finite. 

It is also worth mentioning that the localization length 
scales, approximately, with the mean impurity-impurity 
distance. As Fig. [4] shows, this is clear at the beginning of 
the second and third plateaus, where a linear increase of 
the localization length with energy is found. This linear 
behavior agrees with a prediction discussed by Thouless 
in Ref. [17] . In this reference the relation between localiza- 
tion length and density of states is stressed, showing that 
for a particular density of states, a Gaussian white noise 
model [21], one finds £ ~ E. The discontinuities in local- 
ization length (Fig. [4j) could then indicate the existence 
of discontinuities in the density of states of the disordered 
system. Such discontinuities would not be surprising since 
the density of states of the clean quasi- ID system diverges 
as (E — £ m ) -1 / 2 at the beginning of each plateau, when 
E ~ £ m . 

In Fig. [5] we show the probability distributions of the 
conductance G and Lyapunov exponent 7 for a fixed en- 
ergy E = 1.6ft 2 /(m*p 2 ) and three different wire lengths, 
corresponding to 200, 1000 and 5000 impurities and d\ = 
2.5p. Comparison with the prediction of the theory of dis- 
ordered conductors, Eqs. ([3|) and ([4]), shows the deviation 
in the intermediate regime, N = 1000, when L « 1.5£ and 
the wire is neither in localized nor in metallic regime. In 
the limiting cases, however, either the conductance (metal- 
lic) or the Lyapunov exponent (localized) closely follows 
the expected behavior. 

Figures [3] and [4] also show the semiclassical result in 
each case (continuous line). As mentioned in the Introduc- 
tion, the semiclassical approach neglects interference ef- 
fects by directly composing the transmission probabilities, 
as opposed to the phase- coherent composition of transmis- 
sion amplitudes in fully quantum mechanical calculations. 
Notice that if all phases of Eqs. (fT2j) . (pj) . (fT5l) and (fT6l) 
are neglected the formalism does not depend on the im- 
purity positions (2^, 0$), only on its total number J\f. This 
causes Ohmic behavior, characterized by parameter a in 
Eq. (|6|). It is also possible to relate a with the transmis- 
sion probability of a single impurity Indeed, using 
the composition rule of incoherent scatterers we find the 
total transmission T from 



1 1 =A rf 1 1 
T M M 



leading to the explicit expression 

1 1 \ 1 



1 

M 



do + 0.5di 



(18) 



(19) 
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As seen in Fig.HJ a 1 is close to the localization length 
or, equivalently, a is close to the Lyapunov exponent. 



Fig. 5. Probability distributions of the conductance (left pan- 
els) and Lyapunov exponents (right). Ordered by rows the re- 
sults correspond to increasing wire lengths: 200 impurities in 
(a) and (b), 1000 in (c) and (d) and 5000 in (e) and (f). The 
dashed lines show the normal distributions, Eqs. J3} and J4|, 
from the theory of disordered conductors using an arbitrary 
vertical scale. As expected, Gaussian distributions are found 
in the conductance of the metallic regime (a), and in the Lya- 
punov exponent of the localized regime (f). 



However, both values do not coincide in general, making 
the comparison between quantum and classical conduc- 
tance, Eq. (|7j), rather involved. Another information we 
can extract from the semiclassical result regards the dis- 
continuities in localization length as a function of energy. 
Besides Thouless' argument on the density of states given 
above, we also expect similar drops in the semiclassical 
length a -1 from Eq. ([T9]) . Indeed, the dependence of 
on energy is a staircase with rounded edges. This smooth 
energy dependence is broken by the explicit appearance 
of M, the number of modes, in Eq. (fl9|). A sharp increase 
in a occurs when the energy allows a new mode but 
has not yet increased. 



4 Conclusions 

We have studied the quantum transport properties of tubu- 
lar shaped nanowires with randomly distributed impuri- 
ties. For a given wire length our calculations predict a 
sawtooth-like behavior of the conductance as a function of 
energy; conspicuous drops appearing when the energy is 
just barely enough to allow the propagation of a new trans- 
verse mode. The analysis of our results within the scaling 
theory of disordered conductors explains the observed be- 
havior as a general tendency of the wire to evolve from 
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localized towards metallic regimes with increasing energy. 
Abrupt changes in the density of states occur when a new 
mode becomes propagating, leading the wire back towards 
a more localized regime. 

We have also calculated the localization length of the 
disordered wire for energies covering from one to five prop- 
agating modes. Qualitative interpretation in terms of the 
density of states explains the linear-with-energy depen- 
dence at the beginning of the conducting plateaus, as well 
as the discontinuous drops already mentioned. The semi- 
classical description has been also discussed and shown to 
provide a complementary interpretation of the discontinu- 
ous drops. In the metallic regime, the comparison between 
the average quantum conductance and the semiclassical 
one is complicated when the Lyapunov exponent and the 
semiclassical coefficient for Ohmic behavior do not coin- 
cide. 
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A Single impurity scattering 



The formalism of Sec. l2.2l requires the scattering matrix of 
a single impurity as an input. We find this by solving the 
Schrodinger equation for the corresponding open bound- 
ary problem using the quantum-transmitting boundary al- 
gorithm [22] . The wave function fulfills the equation 



h 2 d 2 



h 2 d 2 



2m* dz 2 2m* p 2 dO 2 



V(z,0)j i/>(z,0) = 

Etl>(z,0) , (20) 



where m* is the electron's effective mass and V(z, 9) is 
given by Eq. (|5]) with the impurity at (z,i, 0i) = (0, 0). Far 
from the impurity, in the parts of the tube that play the 
role of contact leads (or contacts for short), the poten- 
tial vanishes and the Hamiltonian separates in longitudi- 
nal and angular (transverse) contributions. The angular 
problem 

h 2 d 2 cj) n (6) 



2m*p 2 dO 2 



yields the set of transverse modes {4> n ->£n}-> 



K{0) = 



2ir 
2m* p 2 



(A„ =0,±1,...) 



(21) 

(22) 
(23) 



Equation ([20]) can be solved using finite differences in 
a Cartesian grid. However, increasing the number of grid 
points becomes very costly and this method lacks accu- 
racy in the angular integrations with oscillating functions 



like those of Eq. ([22]) . Therefore, we have used a mixed 
approach, in which the z coordinate is discretized in a 
uniform grid while the angular dependence is described 
by expanding in the angular eigenfunctions, i.e., 



(24) 



The unknown band amplitudes ip n (z) fulfill coupled-channel 
equations 



h 2 



2m" 
where 



l/>n(z) + (£n ~ E) l/j n (z) = -^2 V nn> (z)^n' (z) 



V nn '{z) 



2ttJ 



dOV(z,6)e 



i(\ n -\ n ,)6 



(25) 



(26) 



In the contacts the band amplitudes take a similar 
form to Eq. ([TO]) . 

M*) = a c , n e is * k ^ z -^ + b c , n e- is ^ z - z ^ . (27) 

This expression is for a propagating channel, for which 
e n < E and k n = ^/2m*(E — e n )/% is a real number. It 
also applies to evanescent ones, s n > E, if we assume in 
this case a Cy7l = and a purely imaginary wavenumber 
k n = i^/2m*(s n — E)/h. Notice that the output ampli- 
tudes can be obtained from the wave function right at the 
contact position, 

b c ,n = ^n(z c ) - a c , n . (28) 

Substituting Eq. ((28]) in Eq. ([27]) we obtain 

2ia CjU sin (s c k n (z - z c )) , (29) 

which is the quantum-transmitting-boundary equation for 
the contacts. 

Equations ([25]) and ([29]) . for the central and contact 
regions, respectively, form a closed set which does not in- 
voke the wave function at any external point. Of course, 
this is not true for any of these two subsets separately, 
since central and contact regions are connected through 
the derivative in Eq. ([25]) and of ^ n (z c ) in Eq. ([29]) . In 
practice, we use a uniform grid in z with n-point formulae 
for the derivatives (n w 5 — 11) and truncate the expansion 
in transverse bands, Eq. ([24]) . to include typically 50-100 
terms. The resulting sparse linear problem is then solved 
as in Sec. O using routine ME48 [T9j . 
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